Superconducting junctions from non-superconducting doped CuC>2 layers 



V. M. Loktev 1 and Yu. G. Pogorelov 2 
1 N.N. Bogolyubov Institute for Theoretical Physics, Metrologichna 14b, Kiev 0,31,34, Ukraine, 
2 IFIMUP, Umversidade do Porto, R. Campo Alegre, 687, Porto 4169-007 

The theoretical approach proposed recently for description of redistribution of electronic charge 
in multilayered selectively doped systems is modified for a system with finite number of layers. 
A special attention is payed to the case of a finite heterostructure made of copper-oxide layers 
which are all non-superconducting (including non-conducting) because of doping levels being beyond 
the well-known characteristic interval for superconductivity. Specific finite structures and doping 
configurations are proposed to obtain atomically thin superconducting heterojunctions of different 
compositions. 
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An interesting area in nanoengineering of materials 
was opened in a series of experiments by Bozovic et al 
[1] on atomically perfect stacks of selectively doped per- 
ovskite layers. These and some other papers [2, 3] mainly 
used periodic multilayered structures where essential new 
electronic effects, as interface SC between nominally non- 
SC layers [3], appeared due to charge redistribution be- 
tween layers and related shifts of in-plane energy bands. 
The basic condition for SC to appear within few per- 
ovskite layers or even in a single layer is that the lo- 
cal density of hole charge carriers occurs within a def- 
inite, rather narrow, interval: p m in > P > Pmax with 
Pmin ~ 0.07 and Pmax ~ 0.2 (carriers per site). The re- 
quired density distribution results from the correspond- 
ing shifts of in-plane energy bands by local Coulomb po- 
tentials. A simple theoretical model for such processes 
was proposed [4], combining a discrete version of Pois- 
son equation for potential with a band-structure modi- 
fied self-consistent Thomas-Fermi charge density. This 
approach gives exact solutions for infinite periodical and 
some other unbounded systems. However recent studies 
[5-7] showed that pronounced modification of electronic 
ground state and related SC transitions can be obtained 
either in stacks of finite (and small) number of layers 
which is quite promising for practical applications in na- 
noengineered composite devices. The following consider- 
ation aims on an extension of the previous model on an 
arbitrary layered system and establishing criteria for its 
optimum SC performance. This line of study can be seen 
as a practical realization of long envisaged Ginzburg's 
program for ultrathin superconducting states [8] . 

Following the same assumptions as in Ref. [4], we ex- 
press local charge density in j-th layer: pj — e (pj — Xj), 
through the densities pj of mobile holes and Xj of ionized 
dopants (e is the elementary charge) and then present the 
potential difference (fj+i —<Pj between neighbor layers as: 

e JJ \l=j+l 1=1 ) 

This value is obtained considering the electric field Ejj+i 
in the j,j + 1 spacer as the geometric sum of fields Ei 
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FIG. f : Schematic of nanostructured system with selectively 
introduced dopants (white circles) into each of N layers of 
La2Cu04- 

emitted by each i-th charged layer: E\ — 2npi/ (e e sa 2 ), 
on both sides of this spacer. Eq. 1 involves the in- 
plane and c-axis lattice constants a and c, and e c s is 
the (static) dielectric constant that effectively reduces the 
Coulomb field in the c-direction. Eq. I would be exact 
for a stack of mathematical planes, with uniform in-plane 
charge densities pj and separation c, and it should be a 
reasonable model for real La2- a; Sr 2 ;Cu04 layers where pj 
delocalized holes and Xj localized dopants are distributed 
in different atomic planes within the period c of jth layer. 
The adopted form of purely dielectric screening is justi- 
fied in neglect of c-hopping processes, accordingly to their 
above mentioned weakness, also this model neglects po- 
larization effects from the insulating substrate. We note 
that the charge densities pj naturally vanish in uniformly 
doped (pj = Xj), including undoped (pj — Xj = 0), sys- 
tems. Otherwise, the hole carrier density pj can be re- 
lated to the local potential ifj using the respective density 
of states (DOS) gj(e): 

fW/2-etpn 

Pi =2 9j(e)de (2) 

J £f 

with the spin factor 2 (this zero-temperature formula is 
justified for all the considered temperatures T < T c ). 
Thus the role of c-hopping in this model is reduced to 
establishing the common Fermi level £f for all the layers. 
Using the simplest approximation of rectangular DOS: 
gj(e) = 1/W within the bandwidth W, we arrive at the 
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linear relation between pj and tpj : 
1- Pj 



W -e F . 



(3) 



Then, inserting Eq. 3 into Eq. 1, we obtain the linear 
relation between carrier and dopant densities, referred to 
the j,j + 1 spacer: 



Pj+i - Pj = 



1=1 



N 
1=3 + 1 
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, (4) 



where the dimensionless quantity: 
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is the single material parameter of the model. Finally, 
subtracting the relations, Eq. 4, for j,j + 1 and j — 
l,j spacers leads to simple linear equations for the hole 
carrier densities in neighbor layers only: 



Pj+i +Pj-i - (2 + a)pj = 



(6) 



The advantage of Eq. 6 against possible analogous rela- 
tions between the potentials (fj is in avoiding the need to 
know the Fermi level position (doping dependent). For 
an infinite stack of layers, summing these equations in all 
j would automatically assure the total electroneutrality 
condition pj — 0, and this was just the way used in 
Ref. [4] to obtain a more detailed alternative to the phe- 
nomenological Thomas- Fermi treatment [6]. However, 
for a finite stack of j = l,...,N layers, this condition 
should be additionally imposed besides the N — 1 rela- 
tions, Eq. 4, in order to completely define all the N 
densities pj. Since Eq. 6 in this case only applies for 
the internal layers j = 2, . . . , N — 1, one needs two more 
equations which can be obtained from Eq. 4 for terminal 
layers, j = 1 and j = N — 1, under the electroneutrality 
condition: 



(1 + —p 2 = ax ± 
-Pn-1 + (1 + o)Pn = ctx N . 



(7) 



Thus the non-uniform linear system of N Eqs. 6, 7 can 
be presented in the matrix form as: 



1 



(8) 



where the finite TV-stack of layers generates the "Lapla- 
cian" matrix: 



L 



/ 1 -1 ... \ 
-12-10 ... 
0-12-1 ... 

... 0-12-1 
V ... 0-11/ 
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FIG. 2: Modulated electronic configurations by shifted en- 
ergy bands for the samples of selectively doped layered sys- 
tems, a) For a finite stack of iV = 6 layers with doping levels 
xi — X2 — %3 = 0.45, Xi = £5 = x& = (light columns) and 
the localization parameter a = 1, the carrier densities (dark 
columns) are calculated from Eqs. 8, 9. b) For an unlimited 
system with periodic repetition of the same stack, those are 
calculated from Eqs. 8, 10. The dashed lines mark the val- 
ues Pmin and Pmax , delimiting the interval of carrier densities 
where superconductivity should exist. 



and the TV-vectors: 



f 




and 



It is seen from Eq. 8 that a plays the role of localization 
parameter: the carrier density strictly coincides with the 
doping distribution in the limit a — > 00, otherwise it is 
spread beyond this distribution, the stronger the smaller 
a is. Generally, the standard solution p = R(a)~x* with 

the resolvent R = (l + a~ 1 Lj gives the densities pj 

in terms of the doping levels Xj and of the localization 
parameter a as for the instance in Fig. 2 with TV = 6, 
a = 1 (a reliable estimate for real La2Cu04 [9? , 10]) and 
xi = X2 — X3 = 0.45, X4 = X5 = £6 — 0. It is of interest 
to compare this solution to that for an unlimited system 
with the same distribution of dopants Xj but periodically 
repeated, so that the L matrix is replaced by its periodic 
version: 
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\ -1 ... -1 2 / 
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A notable difference in the resulting distributions of den- 
sities pj is seen in Fig. 2a, b. An important feature of 
superconducting layers formed in each of these structures 
(4th layer in Fig. 2a and 4th and 6th layers in Fig. 2b) 
is that they are realized in nominally undoped La2Cu04 
and thus can be expected almost free of undesirable disor- 
der effects (as scattering by defects and consequent fluc- 
tuations of the SC order parameter [11]). 

Also a comparison of such calculation with the exper- 
iment data [7] is presented in Fig. 3, using the real val- 
ues of structure parameters: N = 12 and x\ = x 2 = 
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FIG. 3: Comparison of calculated (dark columns) and exper- 
imentally determined (white circles) carrier densities for the 
real system with N = 12, the doping levels xi = xi — X4 — 
x$ — Xfj = 0.45, X7 = 0.17, x$ = 0.06, xg — 0.02, xw = am — 
£Ci2 = (light columns) and the localization parameter a = 1. 



X4 = xs — xq = 0.45, X'j = 0.17, x$ = 0.06, xg = 0.02 
(the three last values are due to the interdiffusion of 
Sr dopants into nominally undopcd La2CuC>4 layers), 
xio = xn = x\2 = and the same localization pa- 
rameter a — 1 as before. The resulting distribution 
of carrier densities (black columns) displays an excel- 
lent coincidence with its measured values. In partic- 
ular, the carrier density occurs within the SC interval 
just in the N = 8 layer with the value p$ « 0.106 that 
just corresponds to the observed transition temperature 
T c « 32 K when used in the phenomenological formula 
[4] T c (p) = (p- p mm ) ( Pmax - p) T* with T* a 9000 K. 

From such a good agreement, one expects that this ap- 
proach can be also used for design of new configurations 
with tailored superconducting performance. In particu- 
lar, an interesting possibility is to build ultra-quantum 
hetero junctions of two types: superconductor- insulator- 
superconductor (SIS) and superconductor-normal metal- 
superconductor (SNS), each component being restricted 
to a single cuprate layer. Such junctions could realize 
a single layer limit of already discussed thicker SIS and 
SNS nanostructures with giant proximity effect [15]. It 
should be noted that since the localization parameter a 
value (see Eq. 6) is rather fixed by the choice of the 
building material (with a « 1 for La2Cu04), the practi- 
cal control parameters in this process must be the total 
number N of layers in the stack and particular doping 
levels Xj in each layer. 



Thus, one possible simple structure to produce a SNS 
junction can consist of N = 5 cuprate layers with 
the doping levels defined for instance as: x\ — x§ = 
0.45, X2 — X3 — X4 = (that is, nominally all non- 
superconducting). From Eq. 8, the resulting carrier den- 
sity distribution: p± — p§ w 0.29, correspond to nomal 
(overdoped) metal layers, p2 — Pi ~ 0.12 to supercon- 
ducting layers with a high enough transition tempera- 
ture T^ lgh w 37 K, separated by the layer with p 3 « 0.08 
and low transition temperature T l c ow m 11 K. Then, in 
the temperature range T l ° w < T < T^ lgh one should 
obtain a SNS heterostructure. If so, the quasiparticle 
spectrum of this junction will present a peculiar combi- 
nation of gapped (a kind of size quantization) and gap- 
less branches with interesting IR absorption and electric 
transport properties. 

As to the SIS heterostructure, it is rather difficult to 
be obtained in the suggested N — 5 stack, but it can be 
achieved, e.g., by adding one more undoped layer (i.e., 
from 3 to 4) to the above structure, or introducing elec- 
tronic doping in the central layer (making X3 < 0). Un- 
like the above SNS case, the resulting SIS junction should 
display a quasiparticle spectrum with gapped branches 
only. 

In conclusion, an extension of the recent electrostatic 
model for charge redistribution in non-uniformly doped 
multilayered systems is proposed for finite (mostly small) 
number of layers. Formal solutions of this model are 
mainly analyzed in the parameter range actual for the 
experimentally investigated L^-xSrajCuC^ multilayers. 
The distinctions of finite stacks from previously stud- 
ied unlimited or cyclic systems are indicated. Some 
new specific arrangements of doped and undoped layers 
are suggested for realization of artificial atomically thin 
heterostructures with unusual electronic excitation spec- 
trum, potentially interesting for applications in nanocom- 
puting devices. Such artificial structures may present 
also an interest for their behavior under applied mag- 
netic field. 
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